home
***
CD-ROM
|
disk
|
FTP
|
other
***
search
/
The Datafile PD-CD 1 Issue 2
/
PDCD-1 - Issue 02.iso
/
_utilities
/
utilities
/
001
/
meschach
/
!Meschach
/
c
/
machine
< prev
next >
Wrap
Text File
|
1994-08-03
|
4KB
|
175 lines
/**************************************************************************
**
** Copyright (C) 1993 David E. Stewart & Zbigniew Leyk, all rights reserved.
**
** Meschach Library
**
** This Meschach Library is provided "as is" without any express
** or implied warranty of any kind with respect to this software.
** In particular the authors shall not be liable for any direct,
** indirect, special, incidental or consequential damages arising
** in any way from use of the software.
**
** Everyone is granted permission to copy, modify and redistribute this
** Meschach Library, provided:
** 1. All copies contain this copyright notice.
** 2. All modified copies shall carry a notice stating who
** made the last modification and the date of such modification.
** 3. No charge is made for this software or works derived from it.
** This clause shall not be construed as constraining other software
** distributed on the same medium as this software, nor is a
** distribution fee considered a charge.
**
***************************************************************************/
/*
This file contains basic routines which are used by the functions
in meschach.a etc.
These are the routines that should be modified in order to take
full advantage of specialised architectures (pipelining, vector
processors etc).
*/
static char *rcsid = "$Id: machine.c,v 1.4 1994/01/13 05:28:56 des Exp $";
/**** For Acorn Desktop C ****/
/* Naughty poking about to get an isatty function */
/* I don't guarantee that this works */
#define __system_io
#define FILEHANDLE unsigned
#include <stdio.h>
#include "machine.h"
int isatty(int inf)
{ return ((FILE *)inf)->__file == 0;
}
/**** isascii - or at least my guess at it ****/
#include <ctype.h>
int isascii(char *s)
{ for (;*s != '\0'; s++)
{ if (!isprint(*s))
{ return 0 /*FALSE*/;
}
}
return 1 /*TRUE*/;
}
/* __ip__ -- inner product */
double __ip__(dp1,dp2,len)
register Real *dp1, *dp2;
int len;
{
#ifdef VUNROLL
register int len4;
register Real sum1, sum2, sum3;
#endif
register int i;
register Real sum;
sum = 0.0;
#ifdef VUNROLL
sum1 = sum2 = sum3 = 0.0;
len4 = len / 4;
len = len % 4;
for ( i = 0; i < len4; i++ )
{
sum += dp1[4*i]*dp2[4*i];
sum1 += dp1[4*i+1]*dp2[4*i+1];
sum2 += dp1[4*i+2]*dp2[4*i+2];
sum3 += dp1[4*i+3]*dp2[4*i+3];
}
sum += sum1 + sum2 + sum3;
dp1 += 4*len4; dp2 += 4*len4;
#endif
for ( i = 0; i < len; i++ )
sum += dp1[i]*dp2[i];
return sum;
}
/* __mltadd__ -- scalar multiply and add c.f. v_mltadd() */
void __mltadd__(dp1,dp2,s,len)
register Real *dp1, *dp2;
register double s;
register int len;
{
register int i;
#ifdef VUNROLL
register int len4;
len4 = len / 4;
len = len % 4;
for ( i = 0; i < len4; i++ )
{
dp1[4*i] += s*dp2[4*i];
dp1[4*i+1] += s*dp2[4*i+1];
dp1[4*i+2] += s*dp2[4*i+2];
dp1[4*i+3] += s*dp2[4*i+3];
}
dp1 += 4*len4; dp2 += 4*len4;
#endif
for ( i = 0; i < len; i++ )
dp1[i] += s*dp2[i];
}
/* __smlt__ scalar multiply array c.f. sv_mlt() */
void __smlt__(dp,s,out,len)
register Real *dp, *out;
register double s;
register int len;
{
register int i;
for ( i = 0; i < len; i++ )
out[i] = s*dp[i];
}
/* __add__ -- add arrays c.f. v_add() */
void __add__(dp1,dp2,out,len)
register Real *dp1, *dp2, *out;
register int len;
{
register int i;
for ( i = 0; i < len; i++ )
out[i] = dp1[i] + dp2[i];
}
/* __sub__ -- subtract arrays c.f. v_sub() */
void __sub__(dp1,dp2,out,len)
register Real *dp1, *dp2, *out;
register int len;
{
register int i;
for ( i = 0; i < len; i++ )
out[i] = dp1[i] - dp2[i];
}
/* __zero__ -- zeros an array of floating point numbers */
void __zero__(dp,len)
register Real *dp;
register int len;
{
#ifdef CHAR0ISDBL0
/* if a floating point zero is equivalent to a string of nulls */
MEM_ZERO((char *)dp,len*sizeof(Real));
#else
/* else, need to zero the array entry by entry */
int i;
for ( i = 0; i < len; i++ )
dp[i] = 0.0;
#endif
}